CRITICAL BRANCHING RANDOM WALKS WITH SMALL DRIFT 



XINGHUA ZHENG 

Abstract. We study critical branching random walks (BRWs) 17 W on Z + where for 
each n, the displacement of an offspring from its parent has drift 2/3/y / n towards the origin 
and reflection at the origin. We prove that for any a > 1, conditional on survival to genera- 
tion [n a ], the maximal displacement is asymptotically equivalent to (a — l)/(4/3)y/rilog n. 
We further show that for a sequence of critical BRWs with such displacement distributions, 
if the number of initial particles grows like yn a for some y > and a > 1, and the par- 
ticles are concentrated in [0, 0(y/n)], then the measure-valued processes associated with 
the BRWs, under suitable scaling converge to a measure-valued process, which, at any 
time t > 0, distributes its mass over M. + like an exponential distribution. 



1. Introduction 

Durrett et al. ( 1991 ) and iKesten (1995) studied the maximal displacement of critical 



branching random walks (BRWs) on the real line conditioned to survive for a large number 
of generation s. When the spatial displacement distribution has drift \i > 0, the results in 
Durrett et al.l (|199l|) imply that conditional on the event that the BRW survives for n gener- 
ations, the maximal displacement of a particle from the position of the initial particle will 
be of order Op(n). The main result in Kestenl ( 1995 ) asserts that if the spatial displacement 
distribution has mean and finite (4 + 5)th moment, then conditional on the event that 
the BRW survives for n generations, the maximal displacement will be of order Op(^Jn). 
The sharp difference between these two results gives rise to the following natural question: 
What happens if the spatial motions have "small drift"? 

In this paper we supplement these results by showing what happens for BRWs on the 
nonnegative integers Z + with small negative drift and reflection at 0. Assume that is 
a sequence of critical BRWs on the half line Z + = {x E 7h : x > 0}, each started by 
one particle at the origin, that evolve as follows: (A) At each time t = 1,2, ... , particles 
produce offspring particles as in a standard Galton- Watson process with a mean 1, finite 
variance cr 2 offspring distribution Q. (B) Each offspring particle then moves from the 
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location of its parent according to the transition probabilities P = P^> n ), where /3 > 0, 
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for x > 1; 



for x > 1; 



The spatial motion is hence slightly biased towards the origin, which serves as a reflecting 
barrier. Such a BRW can be used to model, for example, a branching process occurring 
in a V-shaped valley, where the particles, due to gravity, have a slight tendency to move 
towards the bottom. In lKac ( 1947|) the afore-described slightly biased random walk is used 
to model the mo tion of "he avy Brownian particles" in a container with its bottom as a 
reflecting barrier. Kac ( 1947h also states about the reflecting barrier that "the elucidation of 
its influence on the Brownian motion is of considerable theoretical interest". In this article 
we will study the influence of the barrier on the BRW. 

Denote by (x) the number of particles in the nth BRW at location x at time 
t, and by R\ ' the location of the rightmost particle at time t. Our main interest is in the 
conditional distribution of given that the process survives for [n a ] generations. 
For a < 1, the effect of the drift —2f3/-Jn will be negligible compared to diffusion effects 
over this time interval, and for a = 1 it is just large enough to match the diffusion effects. 
Thus, we will focus on the case when a > 1. 

Theorem 1. When > 0, for each a > 1 and e > 0, the range -R^li at time [n a ] satisfies 
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where for any k G Z + , 

(3) = {U^ survives to generation k}. 

It is natural to consider in connection with the behavior of the maximal displacement 
the process-level scaling behavior of the BRWs. To this end, consider a series of BRWs 
{X^} on the set Z + of nonnegative integers that evolve by the rules described above, but 
with arbitrary initial states . (In Theorem[T]the initial state consisted of a single particle 
located at the origin 0.) For integers x, k > 0, set 



(4) 



# particles at x at time k. 



For any subset / C R + , let 



Finally, let 
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Recall that by Kolmogorov's estimate for critical Galton- Watson processes (see (|2"TT) be- 
low), if the nth BRW X^ n > is initiated by 0(n a ) particles, then the total lifetime of the 
process will be on the order of Op(n a ) generations. If a < 1, then the effect of the drift 
over a time interval [0, 0(n a )] is too small to be felt. If a = 1 then the drift will be just 
large enough to be felt, and so for large n the BRW X^ n \ suitably rescaled, will look like 
a Dawson- Watanabe process on the halfline [0, oo) with drift —2(3 and r eflection at (fo r 
th e convergence of o rdinary BRWs to Dawson-Watanabe processes, see IWatanabd (| 1968). 



or 



Etheridge (2000); Perkins The case we will focus on is again when a > 1, as 



in this case the effect of the reflecting barrier at dominates the diffusion effects over the 
lifetime of the branching process, and the result is an entirely different scaling behavior: 

Theorem 2. When (3 > 0, assume that for some a > 1, 



,(n) 



(5) 



— > y > 0, as n — > oo, 



and {Xq (yn-)/n a } n >i is tight, i.e., for any e > there exists C > such that for all n, 



(6) 



< e. 



-(n) 



Then the measure-valued processes \X^L t \{\/n-) / 'n a : t > 0J converge, in the sense of 

convergence of finite-dimensional distributions, to a process (X t : t > 0), where (X t ) t > Q is 
such that for allt > and < a < b, 

(7) X t ((a, b)) = Y t ■ (exp(-4/3a) - exp(-4/36)) := Y t ■ vr((a, b)). 
Here Y t is the Feller diffusion: 

(8) dY t = a^/Y t dW t , Y = y. 

Observe that we do not require the initial measures X^ ( y/n-) / n a to converge; what we 
only require are (i) the total mass converges, and (ii) the particles are not too spread out. In 
particular, we cannot guarantee that Xq (^/n-)/n a =^ X . Theorem [2] says that one has 
finite dimensional convergence on (0, oo). 

The Feller diffusion (Y t ) defined by © is the limit of (zfcL/n a ): 



(9) 



y(n) 

n a 



{Y t ) o nJ D([0,oo) 



see lFellen (Il939lll951|) . See Chapter XI of iRevuz and Yon (|1999|) for some basic properties 
of the Feller diffusion. The limiting process X t hence can be described in this way: its 
total mass evolves like the Feller diffusion Y t , but the distribution of th e ma s s Y t a t any 
time t > is always the exponential distribution it. As is proved in iKacl dl947h . the 
exponential distribution n is the stationary distribution of a diffusion process obtained by 
suitably normalizing the RWs defined by (OQ) and taking limit as n — > oo. 
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The following elementary relation between the expected number of particles at a site y 
in generation m for a critical BRW and the m-step transition probability P(S m = y) of the 
random walk will be frequently used: if the critical BRW is started by one particle at site 
x, and U m (y) stands for the number of particles at site y in generation m, then 

(10) EU m (y) = P(S m = y\S = x). 

This is easily proved by induction on m, by conditioning on the first generation and using 
the fact the the offspring distribution has mean 1 . 

The structure of this article is as follows: in Section 2 we prove some properties of the 
random walks on the half line, in Section 3 we prove Theorem [Q Theorem [2] is proved in 
Section 4. 

Notation. We follow the custom of writing / ~ g to mean that the ratio f/g converges 
to 1. For any a, b G R, a A b : = min(a, b) and a V b := max(a, b). Throughout the paper, 
c, C etc. denote generic constants whose values may change from line to line. For any 
x > 0, [x] denotes its integer part, i.e., the greatest integer no greater than x. The notation 
Y n = op(f(n)) means that Y n /f(n) — > in probability; and Y n = P (f{n)) means that 
the sequence \Y n \/ f(ri) is tight. 

2. Random Walks 

Throughout this article we use the notation {S m } m > = {Sm' n ^} to denote a random 
walk with transition probabilities P = p(^' n ) defined by equation (OQ); use {Smj^o to 
denote the simple random walk on Z + with reflection at 0; and use {S m } m >o to denote the 
simple random walk on Z. Furthermore, for any such random walks, e.g., {S m }, for any 
x, y G Z + and m G N, P x {S m = y) = P(S m = y \ S = x) is the probability that S m 
started at x finds its way to site y in m steps. 

The following lemma says that the random walk S m which has drift towards the origin 
is stochastically dominated by the reflected simple random walk S m . 

Lemma 3. For any (3 > 0, n G N and x G Z + , we can build random walks {S m } m > ~ 
p(^' n ) and {S m } m > ~ P on a common probability space so that 

So = Sq = x, and S m < S m , for all m. 

Proof. It suffices to prove the result for the case x > 0; the case x = then follows since 

S 1 = S 1 = 1^ 

Let S = S = x. At time 1 sample a U x ~ Unif (0, 1). If U x < 1/2 + P/y/n, then let 
Si = x — 1, otherwise let Si — x + 1. In the meanwhile, ifUi < 1/2, then let Si — x — 1, 
otherwise let Si = x + 1. Clearly {S Q , Si} and {S 0: Si} follow their laws respectively and 
Si < Si. Now suppose that we have built {S m } and {S m } up to time m, and we have 
S m < S m . If S m < S m , we must have S m < S m — 2 since at each step the difference 
between the jumps is either or 2; now because at each step the random walks can at most 
jump 1, at time m + 1, we must still have S m+ i < S m+ i. In the other case when S m = S m , 
if S m > then we can build S m+ i < S m+ i just as at time 0; otherwise S m = 0, then 
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necessarily S m+ i = S m +i = 1. Thus, we have proved that we can build {S m } and {S m } 
up to time m + 1. By induction, the conclusion holds. □ 

Lemma 4. For any fceN, any x > k, and any m > 0, 

(11) P x (S m >x + k) <P°(max|5<| > k). 

i<m 

Moreover, there exist C > and b > such that 

(12) P° (max. ISA > k] < Cexp ( -—"] , for all m. 

\%<m J \ m J 

Proof. Inequality (fTTj) holds because in order that S m > x + k, either the random walk 
{Si}i< m nas never visited 0, in which case it just evolves like a simple random walk whose 
maximal deviation from x is no less than S m — S > k, or the random walk {Si}i< m 
has visited in which case it evolves like a simple random walk before hitting 0, and the 
maximal deviation from x before time m is no less than x > k. 

Now let us prove (fl"2l) . First recall the fact that for the simple random walk {S m \ S = 
0}, there exists b > such that 



12 



(13) supi?exp b := C < oo, 

m 



see, e.g., Exercise 2.6 in lLawler and Limid (I2007|) . Now by the submartingale maximal 
inequality, we get 

P° > *) = P- (« exp( 9 |S,n > «p(«0) < ^g|ff • 

Inequality ([121) follows by taking to be 6/m and using (fl"3T) . □ 

Next lemma indicates that if two random walks and S 1 ^ have the same drift 2(3 / y/n 
towards the origin, and are such that Sq — Sq is a positive even number, then is stochas- 
tically dominated by S^. 

Lemma 5. For any fixed (3 > 0, n E N, < %\ ^ i 2 , and a random walk {^} m >o ~ P^'™) 
with Sq = 2ii, we can build a coupling random walk {Sm}m>o ~ P^'™) with Sq = 2i 2 on 
a possibly extended probability space such that 

S]n — ^m: f° r a ^ m i tf^ < ^ 

Sjn > S^, for all m, ifi\ > i 2 

Similar conclusion holds if we change the initial positions of {S^} and {S^} to Sq = 
2/; • L,Sg 2i 2 ■ 1. 

Proof. We shall only prove for the case where Sq = 2ii, Sq = 2i 2 and i\ < i 2 . We will 
build {S^} step by step: if S^ > 0, then S^ +1 moves in the same direction away from S^ 
as S} n+1 does, i.e., 

C-2 q2 , (ql _ ql 

w m+l w m ' V w m+1 '-'mil 



(14) 
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otherwise if = 0, then choose S^ +1 according to distribution (QQ). Since Sq — Sq = 
2{i 2 — ii) is even and at each step the difference between the jumps is either or 2, the two 
random walks cannot cross each other and will either never meet, or merge after they meet. 
The dominance (fl4)) follows. □ 

We now loo k more closely at the random walks {S m } ~ P = P^> n ) . Based on the results 



in 



KacNl947h we show the following. 



Proposition 6. For any fixed (3 > 0, a > 0, and any nonnegative integer sequences {s n }, 
{m n } with s n = 0(^Jn) and lim n m n /(n(\ogn) 2 ) > 0, the random walks {Sm^ \ = 
s n } ~ P (/3 ' n) satisfy 

(15) lim P (S<£ > aV^ | 4 n) = s n ) = exp(-4/3a), 
and 

p(stl>a^Elogn\S^^s n ) 

(16) lim — v L = i. 



Proof. When a = 0, (|15l) and (|16l) clearly hold. So below we assume that a > 0. 
Let 

ff = gW = I--^ = , and p = p(») = i + -^. 

2 vn 2 vn 



By (41) inQ dl947h . for any jfe > 0, 



(17) 



A' 



P 1 I 1) (l _|_ C_1^n+fc+mnN 

2pg W v ; 



:=p* (k) + R mn (k), 



f0) = COS ie-2^L s -^, ee[o,7r}. 

\/n sinf 



where for any z > 1, 



We first estimate the main term p* mn {k). Depending on whether s n + m n is even or odd, 

Sml only takes even or odd values. We shall only deal with the case when s n + m n is even. 
In this case, 



h 



k>a^n k>a^/n, k even 

Using the sum formula for geometric series and noting that 

1/3 
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one can easily show that 

(18) lim V p* mn (k)=exp(-4f3a). 
Similarly, 

ayfri log n Prn n 

(k) , 

(19) hm ~ v _ 4flQ = 1. 

It remains to show that the remainder terms R mn (k) decay rapidly as n — > oo. In fact, by 
the simple bound 

| sin i9\ < i sin 0, for all 9 e [0, 7r], 

we get 

1/^)1 <1 + ^- 



Hence, since s n = 0(^/n), 

/ \ k/2 

\R mn (k)\ <C(2y/pq) mn - (J J (1 + 2*0) 

< Cexp(-2/3 2 m n /n) • exp(-kp/y/n)(l + 2Jfe£). 

As 

oo 

J2 e M-k(3/Vn)(l + 2kf3) = 0(y/n), 

k=l 

and lim„ m n /(n(log n) 2 ) > 0, CES) and (QB follow from CD2) and (Tr9l . □ 

3. Proof of Theorem CD 

We first recall some well known facts about critical Galton- Watson processes. Let 
a 2 < oo be the variance of the offspring distribution Q. Then, if Z m is the number of 
particles at time m with Z — 1, and G m = {Z m > 0} is the event that the Galton- Watson 
process survives to generation m, then 

(20) Var(Z m ) = ma 2 , 

(21) Pm : = P(G m ) ~ 

1 (7^ 111 

(22) £(Z m |C7 m ) = — , and 



(23) C 



G m =>■ Exp(cr /2), as m — > oo, 



7 

m 

see, e.g., sections 1.2 and 1.9 of lAthreya and Nevl dl972h . Relation (gfl is known as Kol- 



mogorov's estimate; (1231) is Yaglom's theorem. 

We will decompose the proof of Theorem Q] into two steps. In Proposition|9]we show that 
for any £ > 0, (« — l—£:)v / nlogn/(4/3) is an asymptotic lower bound for i?r™^,. Proposition 
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[10]says that (a — 1 + e)y/n\ogn/ (4/3) is an asymptotic upper bound. Theorem Q] follows 
by combining these two propositions. 

To prove Theorem [TJ we will follow the strategy used to prove Theorems 4 and 5 in 



Lalley and Zhengl (|2007r) . namely, changing the conditional event G k to some event defined 



with respect to a generation m(k) < k. 



The following two lemmas are Lemmas 17 and 18 in Lalley and Zheng (2007) 



Lemma 7. Suppose that on some probability space (Q, J 7 , P) there are two events Ex, E 2 
with P(Ex)P(E 2 ) > such that 

P(E 1 AE 2 ) 

(24) —-— — < e 

where ExAE 2 is the symmetric difference of Ex and E 2 . Then 

(25) \\P(-\Ex)-P(-\E 2 )\\ TV <2e, 

where P{-\Ej) denotes the conditional probability measure given the event E it and \ \ ■ \ \tv 
denotes the total variation distance. 

Lemma 8. Let m(k) < k be integers and e k > be real numbers such that m(k)/k — > 1 
and e k — > as k — > oo. Then 

(26) lim P( p^f fc) = 0, 
where 

G{k) = {Z k > 0} and H{k) = {Z m[k) > ke k }. 

By Lemmas [7] and [H we can change the conditioning event G k = {Z k > 0} to H k = 
{Z m (k) > ke k }, and it suffices to prove the convergence in Theorem \T\ when the condition- 
ing event is H k rather than G k . The advantage of this is that, conditional on the state of 
the BRW at time m(k), the next k — m(k) generations are gotten by running independent 
BRWs for time k — m(k) starting from the locations of the particles in generation m(k). 

We now show that (a — 1 — e) y/n log n / (4/3) is an asymptotic lower bound for R^L . 



Proposition 9. For any e > 0, 

' Gtl I = I . 



lim P ( > ^ -y/nhgn 



Proof. By Lemmas [7] and [8l we can change the conditioning event from G&L to {^™2,_ nL(n) > 



[n a /L(n)}} for L(n) := [(logn) 2 ], where for any k > 0, is the number of particles at 
generation k for the nth BRW . Conditioning on {Z^L nL , n s > [n a / L(n)]}, there will 
be at least X ~ Bin([n Q /L(n)], p n L{n)) number of particles at time [n a ] — nL(n) whose 
families will survive to time [n a ] . For any such particle, among its descendants at time [n a ] 
we uniformly pick one, then the trajectory of the chosen particle from time \n a ] — nL(n) 
to [n a ] will be a random walk following the law F^' n , starting at the location of its an- 
cestor at time [n a ] — nL(n). In this way we get at least Bin([n Q /L(n)], p n L(n)) number 



nL(n) 
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of independent random walks. We would like to show the probability that the maximum 
of the end positions of these random walks is bigger than (a — 1 — e) -^/n log n/ (4/3) is 
asymptotically 1. By Lemma[5l this probability is not increased if we assume that all these 
random walks are started at or I, depending on whether [n a ] — nL(ri) is even or odd. 
But since the random walks have nL(n) steps to go, by relation (fl6l) . no matter whether 
the starting point is or I, for large n, the probability that each random walk is to the 
right of (a — 1 — e)/(4/3) • y/n\ogn is asymptotically n~^ a ~ 1 ~ s \ However we have at 
least X ~ Bin([n°/L(n )1, p n u n )) number of i.i . d. tria ls, and by relation(l2~TT) and Chernoff 
bound (|Chernoffl ( 1952 ) or lAngluin and Valiant! ( 1979h ). for all n sufficiently large, 



/ 1 n a 2 \ / n a 2 1\ 

(27) p \ x -rW)^^)- e ^VW)^^'v)^^ 

It follows that the probability for the maximum of the end positions of these random walks 
to be bigger than (a — 1 — e)/ (A/3) ■ ^fn log n is asymptotically 1 . □ 

Proposition [9] gives the desired lower bound. We now prove the upper bound. 

Proposition 10. For any e > 0, 

2& p (*m 2 '^w 1 ' ^ los " I G ' : "') = 1 

Proof. For any e n ->■ 0, define = {^[ ( "2]_ n > {[n a ] - n) ■ e n }. Applying Lemmas [7] 
and [8] once we see that we can change the conditioning event from Gj™L to H^h; applying 

(n) 

these lemmas again we see that we can change the conditioning event to G^_ n . Since 



a > 1, by relation (fl6l) . the probability that each random walk is to the right of (a 
1 + e/2)/ (4/3) • y/nlogn at time [n a ] — n is asymptotically n~( a ~ 1+e l 2 \ Thus, using 
relations (flOl) and (l22~l) . the conditional expectation of the number of particles to the right 

of (a — 1 + e/2)y/nlogn/ (4/9) in generation [n a ] — n is 

E (Z$ a] _ n | G^_ n ) ■ P ( 5 [nQ ]-n > (« - 1 + e/2)/(4/3) ■ v^logn) 

^(M-^) „-(a-l +£/2 ) ^ 

2 2 
However, by relation (12~TI) . the probability that a Galton-Watson process survives to time n 
is ~ 2/ {na 2 ), hence the number of particles to the right of (a — 1 + e/2)/ (4/3) • y^nlog n 
in generation [n a ] — n whose families survive to time has expectation asymptotically 
equivalent to rr e l 2 , which goes to 0. Therefore if we denote by 

=the rightmost location in generation \n a ] of the descendants of the particles 
which are to the left of (a — 1 + e/2)/ (4/3) • y/n\ogn in generation [n a ] — n, 
then it suffices to show further that 



(28) P > (a - 1 + £)/(4/3) • v^log 



?? 
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By Lemma [51 this probability is not decreased if we assume all the particles to the left of 

(a — 1 + e/2)/(4/3) • y/nlogn at time [n a ] — n are located at M n , where 

the biggest even number < (a — 1 + e/2)/(4/3) ■ -y/nlogn, 

if [n a ] — n is even; 
the biggest odd number < (a — 1 + e/2)/(A/3) ■ y/nlogn, 

if — n is odd. 

In either case, in order that > (a — 1 + e)/ (4/3) ■ y 7 ™ log n, since the ancestors are to 
the left of (a — 1 + e/2)/ (4/3) • y/n log n, at least one descendent will have to travel to the 
right at least s/(8(3) ■ v^log n distance. Hence, since the BRW is critical, we get 

l + e)/(4/3)- v ^logn|G[2 ] _ 



P ( R'tl > (a 



(29) 



< E{z\Xn I G&_J • P M "(S n >M n + e/(8f3) ■ v^logn). 



n 



By Lemma [3l 

(30) P M " (S n >M n + 5/(8/3) • v^logn) < P M ™ (s n > M n + e/(8/3) • v^log 

When n is sufficiently large, M n will be bigger than e/(8/3) ■ \/n log n, so by Lemma|4]we 
get that the probability on the right side of (1301) is bounded by C exp {—be 2 / (64/3 2 ) • (log n) 2 

Using (1291) . noting that E(Z^L_ n \G^l,) = 0(n a ) only grows polynomially in n, we 
get dH. ' □ 



4. Proof of Theorem [2] 

We start with a simple observation. The following lemma about the probabilities of 
survival is a supplement to the convergence in ©. 

Lemma 11. For the total mass processes (Zr"2 t i)t>o and the Feller diffusion (Y t ) t >o, the 
following convergence holds: 

(31) P (zfcl t] > 5n°^j ^P{Y t >5), for all 5 > and for all t > 0. 

Proof. For any t > 0, the convergence in (1311) when S > follows from the marginal 
convergence Z^l t Jn a =>■ F 4 and that P(Yj = 5) = (for any fixed t > 0, by (12TI) and 



(1231 it is easy to show that the marginal distribution of Y t can be described as a Poisson 
sum of exponentials, see, e.g., page 1 36 in Perkins! ( 2002 ). hence is continuous on (0, oo); 
see also page 441 in lRevuz and Yon (|1999|) for an explicit density formula). It remains to 
show 



P(z$ t] >v)^P{Y t >0). 



In fact, by the independence between the BRWs engendered by different initial particles, 



P Z 



r(n) 
[n a t] 







1 



P[n a t], 
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where p rn , as defined in (|2TI) . is the probability that a Galton-Watson process started by a 
single particle survives to generation m. By (TUT) and ©, 

UW^espf-^.z'"')^, ,(-*). 



The right side equals P(Yf = 0), see, e.g., equation (II. 5. 12) in lPerkinsI (|2002h . □ 



Proof of Theorem^ A. Convergence of Marginal distributions. We will show that for 
any fixed t > 0, on the Skorokhod space D([0, oo); M), 

(32) | [nt] [ n V =► (X t ([0, a]) = y ( • tt([0, a])) a>0 . 



a>0 



Let L(n) := [(logn) 2 ], and write 

^[O^a]) _ z{:! t] _ nL(n) xg^O^gaD _^ 

n a n a z (n) ^Kt]-ni W >°}' 

^[n a t-nL(n)] 

For any a > and 5 > 0, we will show the following law of large numbers: 

(33) ^°' fl l) '^U,^^ ' 

\ ^[n»t]-TiL(ri) / 

Claim: If this holds, then we have the finite-dimensional convergence below: for any k G N 

and any < a x < . . . < a k < oo, 

( 34 ) I ^ 1 =^( y *'T([°. a i]))a ll ...,a fc - 

\ / ai,...,aj. 

Note that the LHS and RHS of d32l ) are both increasing p rocesses and the RHS is con- 



tinuous, by Theorem VI. 3. 37 in Jacod and Shiryaev ( 2003|) . the above finite-dimensional 



convergence implies the convergence (1321) as processes on [0, oo). 

We now prove the claim, which is a direct consequence of Lemma ITTl ©, Slutsky's 
theorem and (1331) . We shall only prove the convergence for any single a > 0; the joint 
convergence can be proved similarly. Let / : R — > R be any bounded Lipschitz continuous 
function. We want to show that 

(35) Ef I I ->■ Ef(Y t -n[0,a\). 

In fact, for any e > 0, there exists 5 > such that 

P(0 <Y t <5)<£. 
By Lemma [TTl for all n sufficiently large, 

^(0 < z{:l t] _ nL(n) < Sn°) < 2s. 
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Hence, denote by M = max x \f(x)\, 



< 



+ 



Ef (Xg t] ([0, ^a])/n a ) - Ef(Y t • vr[0, a}) 

/(o) • p (4 n i]-n i( n) = o) - /(0) • p(y t = 0) 



3Me 



J5 / 



7 (n) v (n) 
^[n"f]-nL(n) 



t] — nL(n) 



>8n a } 



-Elf 



'[n a t]-nL(n) 



7r[ °' a] l 1 {<%l-^)>^> 



+ 



/ 1 i {nat] :: L{n) -^m\\ z ^ 



\n a t]-nL(n) 



>5n a } 



-E (f(Y t -n[O ia })-l {Yt>s} 
-I + 3Me + 11 + HI. 



By Lemma [IH I -> 0. By © and Slutsky's theorem, J/J -»■ 0. Finally, /J ->■ by the 
Lipschitz continuity of /, ©, (|33l) and the dominated convergence theorem. 

We now prove the law of large numbers (1331) . by using a mean-variance calculation. Let 

•^Kt]-nL(n) be the configuration of the BRW at time [n a t] - nL(n), Z^l^-riLOn,) be the 

set of particles at time [n a t] — nL(n), and for each particle Ui = uf^ G ^L_ ni ( n )' ^ et 

Xj = a;^ be its location (at time [n a t] — nL(n)), U% l (x) be its number of descendants 
at site x at time k + [n a t] — nL{n), and Z%* be its total number of descendants at time 
k+ [n a t] -nL(n). 

We start with the mean calculation. 



' 7W ' {^t]- n L(n)>*n°} 

^[n-fl-nifn) 



-E 



By relation (fTOl) . 



^(^][Q,^]l4"VnL(n)) i 



^[n«t]-ni(n) 



7 (n) 

P(S nL(n) e [0, A/na] | So = x f ). 



i=i 
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By Lemma [51 if we let 

P°nL(n) P(SnL(n) & [0, \/na\ \ S = 0) 
PnL(n) : = P(SnL(n) £ [0, y/ua] \Sq = \), 



then 

^- { pIim, if^iiseven 

(36) P(S nL{n) e[0,V^a]\S = x t )<{ ™ L(n) . 

PnL(n), if Xi IS Odd. 



Therefore, by Proposition[6]and Lemma ITTL 

E 7(«) ' 1 {<2 t] -„ LM >^} 

(37) \ Z M-nL(n) 



I PnL(n) V PnLin) ' l{z (n) 



>5n a } 



n[0,a] ■ P{Y t > 5). 



On the other hand, by relation (flOl ) again, for any C > 0, 



E ( Z^ X in«t]-nL(n)d C ^ °°)) ) = T£ P ( > I #0 = ^ ) ■ 

' i=l 



= yv H (\n. r^w \ — 1 r> I q, , ... sn.Cic_- ^,( n ) 

nT 

where Zq is the total number of particles at time 0, and x .J is the location of the ith 

initial particle. By the tightness of {Xq % \^-) / n a ] © and (fl"5l) we see for any e > 0, 
there exists C > such that for all n sufficiently large, 

(38) E (^4V^)([^,oc))) < e. 
Therefore by Markov's inequality, 

(39) P (^X$ t] -n H n)(l c ^ oo)) > Vi) < Vi- 

Now by Lemma [5J again, for those particles Ui at time [n Q t] — nL{n) which are to the left 
of Cy/n, if we let 

M n]even = the biggest even number < C\fn\ 
M n - dd = the biggest odd number < C\/n, 
and 

PnL(n) ■= P{SnL(n) £ [0, y/mi] | S = M n , even ) 
PnL(n) ■= P(SnL(n) G [0, ^ | = M„ ;0(W ), 

then 

r . fKT(n)> if ^ is even 

(40) P(S nL(n) E 0, v^a So = O > S J ' . , , 

\PnL(n), if IS Odd. 
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Hence, by Proposition^ Lemma ITTI and (|39l) . 



'E(X$q[0,y/M]\J$ 



■ j l-- nL(n) 

lim inf — 



> lim inf E 



y(n) 



(41) ' 1 {z[:l t] . nL{n) >^; xl:l t] _ nHn) ([C^^))<n^ 

( 7(") _ A? 

-> 1- • f p ( [n"t]-nL(n) V e „ en ocM 

> limint £ p nL(n) A p nL(n) 



> 



A _ £ V ^[0, a ] ■ (P(y t > 5) - y/e) . 



By the arbitrariness of e, we get the desired lower bound 



lim inf E 

n 

>TT[0,a]-P(Y t >5). 



^(^]([0^a])l^ t] -^( w) ) 



^[n a t]-nL(n) 



So, combining it with (1371 ), we get the convergence of expectation 



^[n°t]-n£(n) ; 



ft 



It remains to show that 



^[f»«t]-ni(n) 1 
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By conditioning on ^L nL(n) , we get 



[n a t}~nL(ny 



*r5ti[0,V& 



Var 



n] 



na\ 



-E 



' '[n a t]-nL(n) 
(n) 



7r[0 ' a] i- V%]-^)>^> 



vv 

+ Var 
=/ + II. 



'[n a t]-nL(n) 



\ 



1 {zH 



[n"(]-iiL(ii) 



><5n a } 



- 7r[0, a] ■ 1 { (n) 



>8n a } 



We will show that both terms converge to 0. 

We start with term I. Recall that for each particle m G <2^2 t ]_ n £,( Tl )> U^(x) denotes its 
number of descendants at site x at time k + [n a t] — nL(n), and is its total number of 
descendants at time k + [n a t] — nL{n). By the independence between the BRWs U Ui and 



) 



< 



E 

») 

n a t]- 

E 



(») 

[n a t]-nL(n) 



ii , 6[0,v / na] 



£ [Z Ui T( , 

1 nL(n) 



(n) 



lC [n a t]-nL(n) 



r(n) 



t]-nL(n)( 1 + nL (n)cT 2 ). 



Hence 



I < E 



1 + nL{n)a 2 

y(l%) 

n \n a t\-nL{n) 



->■ 0. 



As to term H, by (1361) . 

^ (^[n^]([°> V^]) I F[n«t]-nL{n)J ^ ^[n^]_ ni ( n ) ' PnL(n) V PnL(n) 

furthermore, on the event {X^ n l t ^_ nL ^ n ~ ) ([Cy/n, oo)) < n a ^}, by <©, 

^ «%]([0, | ^-^»)) ^ ( Z £V„L(») - ^) • PSR) A Klin) 
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Hence, 
II 

'( n ) (\c\ yCT^n I Tr( n ) 



- E l I ^) ^MJ • W- M >^ 



+E (max ((p° L(n) V ^ L(n) - tt[0, a]) 2 , ((1 - ^ A p$ n) - 7r[0, a]) 1 



'V [ %-n i ( n )>^ a ;^ t ]-^(„)([ C7 ^°°))<" a ^} 
=0(Vi), 

where the term y/e in the second inequality comes from <l39l) , and in the last equation we 
used Proposition |6] By the arbitrariness of e, II — Y and hence (l42l) holds. 

B. Convergence of Finite Dimensional Distributions. This follows from the Markov 
property and similar calculations as in Part A. □ 
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